The effect of inter-Landau-band mixing on extended states in quantum Hall system 
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The effect of inter-Landau-band mixing on electron localization in an integer quantum Hall system 
is studied. We find that mixing of localized states with opposite chirality tends to delocalize the 
states. This derealization effect survives the quantum treatment. Extended states form bands 
because of this mixing, as we show through a numerical calculation on a two-channel network 
model. Based on this result, we propose a new phase diagram with a narrow metallic phase that 
separates any neighboring QH phases from each other and also separates each of them from the 
insulating phase. We reanalyzed the data from recent non-scaling experiments, and show that they 
are consistent with our theory. 
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I. INTRODUCTION 

According to the scaling theory of localization 1 ., all 
electrons in a disordered two-dimensional system are lo- 
calized in the absence of a magnetic field. In the pres- 
ence of a strong magnetic field B, a series of disorder- 
broadened Landau bands (LBs) will appear, and an ex- 
tended state resides at the center of each band while 
states at other energies are localized-. The integrally 
quantized Hall plateaus (IQHP) are observed when the 
Fermi level lies in localized states, with the value of the 
Hall conductance, a xy = ne 2 /h, related to the num- 
ber of occupied extended states (n). Many previous 
studiea 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ! 11 : 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 have been focused 
on so-called plateau transitions. The issue there is how 
the Hall conductance jumps from one quantized value to 
another when the Fermi level crosses an extended state. 
There are two competing proposals. One is the global 
phase diagram 4 , based on the levitation of extended states 
conjectured by Khmelnitski: 19 and Laughlin 20 . A crucial 
prediction of this phase diagram is that an integer quan- 
tum Hall effect (IQHE) state n in general can only go into 
another IQHE states »± 1, and that a transition into an 
insulating state is allowed only from the n = 1 state. 
The other is so-called direct transition phase diagram 5 
in which transitions from any IQHE state to the insu- 
lating phase are allowed when the disorder is increased 
at fixed B. So far, most experiments 1 ^ 1 ^ are consistent 
with the direct transition phase diagram although the 
early experiments were interpreted in terms of the global 
phase diagram. 

One important yet overlooked issue regarding IQHE is 
the nature of both plateau-plateau and plateau-insulator 
transitions. In all existing theoretical studies, these tran- 
sitions are assumed to be continuous quantum phase 
transitions. This assumption is mainly due to the early 
scaling experiments 3 . The fingerprint of a continuous 



phase transition is scaling laws around the transition 
point. In the case of IQHE, it means algebraic divergence 
of the slope of the longitudinal resistance in temperature 
T at the transition point. However, recent experiments 18 
showed that such slopes remain finite when the curves are 
extrapolated to T = 0. This implies a non-scaling behav- 
ior around a transition point, contradicting the expecta- 
tion of continuous quantum phase transitions suggested 
by the theories. Thus the nature of these transitions 
should be re-examined. 



The samples used in the non-scaling experiments 1 ^ are 
relatively dirty, and strong disorders should lead to a 
strong inter-Landau-band mixing. In a recent lette»2i, 
we showed that the single extended state at each LB 
center broadens into a narrow band of extended states 
when the interband mixing of opposite chirality is taken 
into account. A narrow metallic phase exists between two 
adjacent IQHE phases and between an IQHE phase and 
the insulating phase. A plateau-plateau transition cor- 
responds to two consecutive quantum phase transitions 
instead of one as suggested by existing theories. In this 
paper we shall present the detailed description of this 
study. 



The paper is organized as follows. The semiclassical 
network model for two coupled LBs is illustrated in Sec. 
II. It is shown that mixing of localized states of oppo- 
site chirality tends to delocalize a state while mixing of 
states of the same chirality does not. Our approach, level- 
statistics technique, is described in Sec. III. In Sec. IV 
our numerical results and discussions are presented, and 
the original data from the non-scaling experiments are 
reanalyzed according to two quantum phase transition 
points in each IQHP-insulator transition. The conclu- 
sions of this paper are summarized in Sec. V. 
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II. THE SEMICLASSICAL MODEL INCLUDING 
INTER-LANDAU-BAND MIXING 

According to the semiclassical theory^, the motion of 
an electron in a strong magnetic field and in a smooth 
random potential can be decomposed into a rapid cy- 
clotron motion and a slow drifting motion of its guiding 
center. The kinetic energy of the cyclotron motion is 
quantized by E n = (n + l/2)fiw c , where lo c is the cy- 
clotron frequency and n the LB index. The trajectory 
of the drifting motion of the guiding center is thus along 
an equipotential contour of value Vo = E — E n , where 
E is the total energy of the electron. The local drifting 
velocity v(r) is determined by (in SI unit) 

v(r) = \jV{r) x B/(eB 2 ) (1) 

where \jV{r) is the local potential gradient. The equipo- 
tential contour consists of many closed loops. Neglect 
quantum tunneling effects, each loop corresponds to tra- 
jectory of one eigenstate. The motion of electrons are 
thus confined around these loops with deviations typi- 
cally order of the cyclotron radius l c = \J h/ (eB). 




FIG. 1: (a) Four neighboring loops in a one-band model for 
the case of Vo < 0. Dashed lines denote quantum tunnelings. 
The arrows indicate the drifting direction, (b) Loops localized 
around a valley and a peak, respectively. The arrows inside 
a loop show the directions of local potential gradient around 
the peak or valley. The arrows on a loop indicate the drifting 
direction. 

To illustrate this semiclassical picture, let us think of 
the smooth random potential as a landscape of many 
peaks and valleys distributed randomly in the plane. 
Imagine that the landscape is filled with water up to 
a surface with the height of value Vo- The equipoten- 
tial contour of value Vo is thus the intersection between 
the land and the water. According to the percolation 
theory^, the percolation threshold of a two-dimensional 
(2D) continuum model is p c — 1/2, where p c is the oc- 
cupation probability of the medium (the land or the wa- 
ter). For simplicity, we suppose that the distribution 
of the random potential is symmetric around zero. By 
symmetry the percolation point of both the land and 
the water is at Vq = in this case. When Vq < 0, 



the occupation probability of land is above 1/2. Thus 
the land percolates and the water forms isolated lakes. 
These lakes are around valleys and their boundaries cor- 
respond to trajectories of localized states. In the case of 
Vo > 0, the water forms a percolating sea and the land 
becomes isolated islands around potential peaks. The 
boundary of each island is an electronic state. In short, 
semiclassical electronic states in a QH system are equipo- 
tential loops. These loops are localized around potential 
peaks for Vo > and around potential valleys for Vo < 0. 
The drifting direction of each loop is unidirectional. This 
means that they are chiral states. From Eq. ^ one can 
see that states around a peak have opposite chirality from 
states around a valley because the directions of the local 
potential gradients around a peak are opposite to those 
around a valley. If one views the plane from the direction 
opposite to the magnetic field, the drifting is clockwise 
around valleys and counter-clockwise around peaks, as 
shown in Fig. ^ Right at Vo = both the land and the 
water percolate, and the intersection between them is the 
trajectory of an extended state. It means that there is 
only one extended state at Vo = for each LB. As Vq 
approaches zero from both sides, the localization length 
£ of the system tends to diverge as 

£<xiK>r (2) 

where the critical exponent v = 4/3 according to the 
classical percolation theory. Quantum effects are ignored 
in the above semiclassical argument. When two spatially 
separated loops on the same equipotential contour come 
close at saddle points of the random potential, quantum 
tunnelings should be considered. An example in the case 
of Vo < is shown in Fig. [Tfa). In the absence of inter- 
band mixing, numerical calculations have suggested that 
there is still only one extended state in each LB while the 
value of the critical exponent v is modified to be around 
7/34 

In the case of strong disorders or weak magnetic field, 
the width of the LBs is comparable with the spacing be- 
tween adjacent LBs (the Landau gap), and inter-Landau- 
band mixings should no longer be ignored. In order to in- 
vestigate the consequences of inter-Landau-band mixing, 
we shall consider a simple system of two adjacent LBs. 
Since we are interested in interband-mixing of opposite 
chirality, we consider those states with energy between 
lower and upper bands which are centered at E\ and E u , 
respectively, as shown in Fig. Ota). Thus, equipotential 
loops are Vi = E - E t > and V u = E — E u < for the 
lower and upper LBs, respectively. Using the semiclas- 
sical theory described in the previous paragraphs, states 
from the upper band should move along equipotential 
loops around potential valleys while those from the lower 
band around potential peaks, as shown in Fig. Hfb). The 
loops for the upper band drift in clockwise direction, and 
those for the lower band in the counter-clockwise direc- 
tion. These two sets of loops are thus spatially separated 
and have opposite chirality. If we assume that the peaks 
and valleys of random potential form two coupled square 
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FIG. 2: (a) Two adjacent Landau bands in the case when 
the disorder broadened band width is comparable with the 
Landau gap. D(E) is the density of states. E u and E\ denote 
the centers of the two bands, (b) Two sets of equipotential 
contours for electronic states of energy E shown in (a) . One is 
from the lower band with Vq = Vi > 0, and the other from the 
upper band with Vo = V u < 0. The ellipses denote the loops 
and arrows on them show drifting directions. The solid curve 
is the schematic plot of the random potential. Two dashed 
horizontal lines indicate two constant potential planes. B is 
the magnetic field. 



lattices, the loops can be arranged as shown in Fig. E^a), 
where P and V denote peaks and valleys, respectively. 
In the absence of interband mixing, the model is reduced 
to two decoupled single-band models and all electronic 
states between the two LBs are localized. If we intro- 
duce interband mixing, the localized loops may become 
less localized. To see that this indeed occurs, let us con- 
sider an extreme case with no tunnelings at saddle points, 
but with such strong interband mixing that an electron 
will move from a loop around a valley to its neighboring 
loop around a peak and vice versa, as shown by B — > C in 
Fig. Ola). Follow an electron starting at A, its trajectory 
will be A — > B — > C — > D —> E ■ ■-. The electron is no 
longer confined on a closed loop, but is now delocalized! 

In the one-band model, an electron can also hop from 
one loop to its neighboring loops by quantum tunnel- 
ings. At a first glance, this effect seems similar to that 
of interband-mixing. However, they are fundamentally 
different. In the one-band model, electronic states for 
a given Vq are of the same chirality. Thus the drifting 
direction of an electron will be inverted when it tun- 
nels into neighboring loops. This means that strong 
tunnelings in a one-band model will induce an effective 
backward-scattering which also localizes the electrons. 
We can understand this by considering a small part of 
the one-band model as shown in Fig. ^a) where all 
loops are moving in clockwise direction. Without tun- 



neling, the trajectory of an electron starting from point 
A is A^B^I^J^ A, a clockwise closed loop. 
With strong tunnelings, the trajectory will tend to be 
A^B^C^D---^H^A, a counter-clockwise 
closed loop. Thus, the tunnelings between loops of the 
same chirality cannot delocalize states. 




FIG. 3: (a) Topological plot of the trajectories of the drifting 
motion of guiding centers (rhombus). The drifting motion 
around a potential peak (valley) is denoted by P (V), and 
their direction are indicated by the arrows. Dashed lines stand 
for interband mixing, and dotted lines for tunneling at saddle 
points. The thick line (A to I) describes the trajectory of an 
electron due to a strong interband mixing, (b) The equivalent 
two-channel network model of (a). Solid and dashed lines on 
each link denote two channels from two LBs. Squares stand 
for saddle points. P, V and arrows have the same meaning as 
those in (a). 

It is worthwhile to explain why we consider only those 
states between two LB centers. For states outside this 
range, both sets of loops are localized around either val- 
leys or peaks. This means that interband mixing mainly 
occurs between two loops localized around the same po- 
sition, and this mixing will not delocalize a state. In fact, 
as explained in the previous paragraph, the mixing of the 
same chirality does not help delocalize an electron. This 
is why we shall consider mixing between spatially sepa- 
rated states with opposite chirality. Of course, it does 
not mean that the mixing of the same chirality has no 
effect at all. As it was found in some previous works 6 , 
this kind of mixing may shift an extended state from its 
LB center. Level shifting due to mixing between states 
of the same chirality may distort the shape of the phase 
diagram, but should not alter its topology. The emer- 
gence of the bands of extended states is exclusively due 
to the mixing between states of opposite chirality. 

Now, we describe our two-channel network model in 
detail. Assume that tunnelings of two neighboring local- 
ized states (loops) of the same band occur around sad- 
dle points, and interband mixing takes place only on the 
links, Fig. |3{a) is topologically equivalent to the model 
shown in Fig. EJb). Fig. OJb) is the schematic illus- 
tration of our two-channel Chalker-Coddington network 
model. It is similar to the model studied in previous 
publications^*^. There are two channels on each link. 
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One, denoted by a solid line, is from the lower LB around 
a potential peak. The other (dashed line) is from the up- 
per LB moving around a potential valley. The arrows in- 
dicate the drifting direction of the two sets of states. At 
each node, the tunneling between two neighboring states 
of the same LB occurs. As shown in Fig. Ufa), let Z^j*. 

and Z 1 ^ be the incoming wave amplitudes of states 1 
and 2 from upper (lower) LB, respectively, and 
and Z^Jn be the outgoing wave amplitudes of the two 
states. The tunneling is described by a SO (4) matrix. 
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where the subscripts u and I denote the upper and 
lower bands, respectively. The elements s^,n and 
are tunneling coefficients of an incoming wave-function 
in the upper (lower) band being scattered into outgo- 
ing channels at its left-hand and right-hand sides, re- 
spectively. s^, n and s^V,i are related to each other as 
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trix. Under quadratic potential barrier approximation, 
— i.e., V(x, y) = —Ux 2 + Uy 2 + V c around a saddle point, 
where U is a constant describing the strength of potential 
fluctuation and V c is the potential barrier at the point, — 
one can show that the left-hand scattering amplitude is 
given bji2& 
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where sin# describes the interband mixing. = 1 ~ 4) 
are random Aharonov-Bohm phases accumulated along 
propagation paths. In our calculations, we shall as- 
sume that they are uniformly distributed in [0, 2n] 2 £. 
In the following discussion, a parameter P, defined as 
v/P/(l + P) = sin 6, is used to characterize the mixing 
strength. P will take the same value for all links in our 
calculations. We hope that this simplification will not 
affect the physics. 
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FIG. 4: (a) A node with four incoming and four outgoing 
channels. Z" 1 ^ is the wavefunction amplitude of the \ th in- 

coming wave from the upper (lower) LB. Z°m'* is that of out- 
going wavefunction amplitude, (b) A link with two channels. 
Z^a) is the incoming (outgoing) wavefunction amplitude 
of the upper (lower) LB. 



where e„ w = [E + V c - + 1/2)E 2 \/E 1 with E 

1 and E 2 = 



the electronic energy, E\ = ^^y/K 



1 with K = 



em 2 



+ 1. The kinetic energies 



of cyclotron motion in the two bands are (n u + l/2)p2 
and (ni + l/2)p2, respectively, where n u m are the band 
indices and An = n u — ni — 1. The dimensionless ra- 
tio E r = E2/E1 = 2 ^/ 1 + approaches 2 from above 

as U or the inverse of ui c increases 2 ^, i.e., the regime of 
strong disorders or weak magnetic field. Since this is the 
regime we are interested in, we choose the value of it to 
be 2.2 in our calculations. For convenience, we choose E 2 
as the energy unit and the cyclotron energy of the lower 
band as the reference point. The energy regime between 
the two band centers is thus E £ [0, 1]. 

Inter-band mixing between two channels on a link as 
shown in Fig. 01b) is described by a U(2) matrix 
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III. THE APPLICATION OF 
LEVEL-STATISTICS TECHNIQUE ON THE 
NETWORK MODEL 

Electron localization length is often obtained from the 
transfer matrix method. For a two-dimensional system, 
however, it is well known that this quantity alone does 
not provide conclusive answers to questions related to 
the metal-insulator transition (MIT)^ 7 .. On the other 
hand, level-statistics analysis has been used in study- 
ing MIT22iS£. Level-statistics analysis is based on ran- 
dom matrix theory (RMT) 3 — . The basic idea is that the 
localization property of an electronic state can be de- 
termined by the statistical distribution function P(s) of 
the spacing s of two neighboring levels. For localized 
states, the distribution is Poissonian Ppe{s) — exp(—s), 
called 'Poissonian ensemble (PE)\ In the case of extended 
states, the nearest neighbor level spacing distribution has 
the following form^S 



P( 



C 1 {(3)s l3 ex P [-C 2 {f3)s' 2 



(7) 
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where C\(j3) and C2(/3) are normalization factors deter- 
mined by J P(s)ds = 1 and J sP(s)ds = 1. The pa- 
rameter (3 is determined by the dynamical symmetry of 
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the system. The case of j3 = 1 is for systems with time- 
reversal symmetry and an integer total angular momen- 
tum and is referred as 'Gaussian orthogonal ensemble'. 
Systems with time-reversal symmetry and a half-integer 
total angular momentum belong to the case of (3 = 4, 
called 'Gaussian symplectic ensemble'. For systems with- 
out time-reversal symmetry (3 = 2, and it is called 'Gaus- 
sian unitary ensemble (GUE)'. 

We shall follow the approach proposed by Klesse and 
Metzler'- 1 . A quantum state of a network model can be 
expressed by a vector whose components are electronic 
wave-function amplitudes on the links. In our case, the 
vector can be written as $ = ({Z^Zf}), where Z l u and 
Z\ are the electron wave-function amplitudes of the up- 
per band (u) and the lower band (I) on the i-th link, 
respectively. As shown by Fertig22, the network model 
can be described by an evolution operator U(E), an in- 
dependent matrix determined by the scattering proper- 
ties of nodes and links in the model. (As an example, 
the evolution operator of a two-channel network of size 
L = 2 with periodic boundaries on both directions is 
constructed explicitly in the Appendix.) In general, the 
eigenvalue equation of the evolution operator is 



U{E)$ a {E) = e^ E ^ a (E), 



(8) 



where a is the eigenstate index of U. The true eigenen- 
ergies {E n } of the system are those energies at which 
Lo a (E) is an integer multiple of 2tt. It has been shown 
by Klesse and Metzler— that the level-spacing statistics 
of the set of quasi- energies {Lo a (E n )} is the same as that 
of {E n }. Thus the localization property of an electronic 
state with an energy E can be obtained by the quasi- 
energies. The advantage of this approach is that all the 
quasi-energies can be used in the analysis so that better 
statistics can be obtained. 

Chalker and CoddingtonSS showed numerically that 
an open boundary condition along one direction creates 
extended edge states along the other direction. In or- 
der to get rid of the edge states, we employ a periodic 
boundary along both directions in our calculation. For a 
two-channel network model of L x L nodes with periodic 
boundaries along both directions, there are 4L 2 compo- 
nents in $. U is thus a (4L 2 ) x (4L 2 ) matrix. However, 
there is a special property of the network model^: the 
nodes scatter electrons only from vertical channels into 
horizontal channels and vice versa. If one separates $ 
into the set of wavefunction amplitudes on the horizon- 
tal links $h and the set of wavefunction amplitudes on 
the vertical links the evolution equation in one time 
step can be written in the following form 



*v(t- 



U V ^h 

Uh^v 6 



*ff(t) 



(9) 



where 6 is the (2L 2 ) x (2L 2 ) zero matrix. Uv^h de- 
scribes how wavefunction on vertical links evolves into 
that on the horizontal links. Similarly, Uh^v describes 
that from horizontal to vertical links. For the detail 



derivation, we refer readers to the example shown in the 
Appendix. The evolution equation in two time steps is 
given as 

$ ff (t + 2) = U v ^ H U H ^v®H{t) (10) 
*v{t + 2) = U H ^vUv^B$v(t) (11) 

Therefore, the evolution matrix in two time steps is 
block-diagonal and the two blocks have essentially the 
same statistical property. We thus need only consider 
one of them. 

We study the model for L=8, 12, 16, 20 and 24. The 
calculation procedure is as follows. Take a realization of 
the random phases, construct the evolution matrix and 
obtain the quasi-energies {u>i}- Put them in descending 
order and calculate the level spacings Sj = (u>i — u>i-i)/5, 
where S is the average of Si. Repeat this procedure for 
sufficient times so that more than 5 x 10 4 level spacings 
are collected for a given E and P. The level-spacing 
distribution function P(s) is thus obtained numerically. 



IV. NUMERICAL RESULTS AND 
DISCUSSIONS 

A. Analysis of the level-spacing spectrum 

In the following, we shall analyze the numerical re- 
sults of the level-spacing distribution function P(s). Our 
purpose is to show evidence for the existence of bands 
of extended states in our model. Due to the chiral na- 
ture of the drifting motion, time-reversal symmetry is 
absent from our semiclassical network model. Then, ac- 
cording to the RMT— , if bands of extended states exist, 
P(s) of them should be the GUE distribution Pgue(s) = 
32ir~ 2 s 2 exp(-4:S 2 /ir). While P(s) of localized states is 
the PE distribution Ppe(s). Since the global shapes of 
GUE and PE are quite different, let us first take a look 
at the global shape of our numerical results of P(s). 
Curves in Fig show P(s) at (E = 0,P = 0.1) (a), 
(E = 0.02, P = 0.1) (b), and (E = 0.5, P = 1.5) (c) 
for L = 8, 12, 16, 20, 24 and comparison with the Gaus- 
sian unitary ensemble distribution Pgue{s), while Fig|B| 
is for (E = 0.0, P = 0.7) (a), (E = 0.02, P = 0.7) (b) 
and (E = 0.5, P = 0.5) (c). The global shape 

of these curves has some common features. All curves 
have a vanishing value when s tends to zero. At small s 
they increase with s and reach a peak at some intermedi- 
ate s. Then they decrease with increasing s and tend to 
vanish at large s. All these features are the same as the 
GUE distribution P GUE (s)2&. Thus most of them are 
close to Pgue(s) at first glance. This raises the ques- 
tion on how to distinguish between extended states and 
localized states by our numerical results. As a simple 
way, it is natural to expect that P(s) of extended states 
should tend closer to Pgue{s) or remain unchanged with 
increasing L while that of localized states should deviate 
from Pgue(s) with increasing L. By careful observation 
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FIG. 5: P(s) vs. s for L = 8, 12, 16, 20, 24 and comparison 
with the Gaussian unitary ensemble Pgue(s). (a) E = and 
P = 0.1; (b) E = 0.02 and P = 0.1; (c) E = 0.5 and P = 1.5. 



FIG. 6: P(s) vs. s for L = 8, 12, 16, 20, 24 and comparison 
with the Gaussian unitary ensemble Pgue{s). (a) E = and 
P = 0.7; (b) P = 0.02 and P = 0.7; (c) E = 0.5 and P = 0.5. 



one can indeed see that curves in each sub-figure of Fig [5] 
tend to be closer to Pgue(s) with increasing L while 
those in FigHJshow the opposite tendency. Thus we can 
use the different tendency of P(s) with increasing L to 
distinguish between extended states and localized states. 
In the following, we shall show quantitatively such op- 
posite tendencies for extended states and localized states 
by considering several characteristic quantities of P(s). 

Let us first consider a characteristic quantity Iq defined 
by Iq = J s 2 P(s)ds/2. It is commonly used to character- 
ize the global shape of P(s) and to exam the localization 
property? 3 . It is well-known that Iq — 1 for localized 
states while Iq < 1 for extended states^. Thus, the fol- 



lowing simple criteria is employed. If Iq of a state with 
energy E increases and approaches 1 with increasing L, 
this state is localized. Otherwise, it is extended. Curves 
in FigElare in vs. mixing strength P at E — (a); 0.02 
(b); and 0.5 (c) for L = 8,12,16,20,24. FigEJa) shows 
that the state at E = 0.02 is localized at zero mixing 
and extended at small P. Then it is localized again after 
P passes a particular P c where Iq of different L cross. 
For the state at the lower band center E = shown in 
Fig|7Jb), it is extended at zero mixing. Then, it shows 
the same feature as the state of E = 0.02 at small and 
large P. Fig[7fc) shows that state at E = 0.5 is always 
localized at small P and extended only for large P(> 1) 
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FIG. 7: Jo vs. P for L = 8,12,16,20,24, (a) £ = 0; (b) FIG. 8: I P (s) vs. s for L = 8,12,16,20,24 and comparison 
£ = 0.02; (c) E = 0.5. with that of P G ue{s). (a) E = and P = 0.1; (b) £ = 0.02 

and P = 0.1; (c) E = 0.5 and P = 1.5. 



where all curves of different system sizes tend to merge 
together. 

It is well-known that a fundamental difference between 
P(s) of localized and extended states is its behavior at 
small s. When s tends to zero, P(s) tends to zero for 
extended states due to level-repulsion while for localized 
states it tends to one due to level- aggregation^. Thus 
we need to consider the behavior of P(s) at small s for 
further test of the results in the last paragraph. It is con- 
venient to consider a function of integrated level-spacing 
distribution at small s defined by Ip(s) = J Q S P(s')ds'. 
The meaning of Ip(s) is the integrated fraction of level- 
spacings smaller than s. Although P(s) in most cases of 
our numerical results is close to the GUE distribution, 



level-repulsion of extended states and level-aggregation 
of localized states should still be expected at small s. 
This leads to the following criteria for localization prop- 
erty : Ip(s) at small s should increase with increasing L 
for localized states while decrease or remain unchanged 
with increasing L for extended states. Thus the behav- 
ior of Ip(s) at small s can serve as another method to 
distinguish between extended and localized states. Fig|H| 
shows Ip(s) for (E = 0,P = 0.1) (a), (E = 0.02, P = 0.1) 
(b) and (E = 0.5, P = 1.5) (c) for L = 8,12,16,20,24 
and comparison with Ip(s) of Pgue(s)- Fig|§| is for 
(E = 0,P = 0.7) (a), (E = 0.02, P = 0.7) (b) and 
(E = 0.5, P = 0.5) (c). One can see clearly that states in 
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FIG. 9: I P (s) vs. s for L = 8,12,16,20,24 and comparison FIG. 10: I P (s = 0.5, P) vs. P for L = 8,12,16,20,24, (a) 
with that of P G ue{s), (a) E = and P = 0.7; (b) £ = 0.02 E = 0; (b) E = 0.02; (c) P = 0.5. 
and P = 0.7; (c) E = 0.5 and P = 0.5. 



FigJSJ show the feature of extended states while states in 
FigOare localized. In order to exam each electronic state 
of fixed electronic energy E in the whole range of mixing, 
we consider Ip(s — 0.5), the fraction of the level-spacings 
less than 0.5. We plot the results of Ip(s = 0.5) vs. P 
at E = 0,0.02 and 0.5 for L=8, 12, 16, 20, 24 in FiglTUl 
Similar with the criteria for Ip(s), we use the following 
criteria. If 7p(0.5) of a state increases with increasing L, 
the state is localized. Otherwise, they are extended. Ac- 
cording to this criteria, curves in Fig^|lead to essentially 
the same results for localization property as obtained by 
the analysis of 7o m Fig[7| consistent with the results in 
the last paragraph. 



Let us now turn to the region of large s. Since Pgu e (s) 
decays faster than Ppe(s) at large s, the behavior of 
P{s) in this region can also be used to show differ- 
ence between extended states and localized states. At 
this region it is convenient to consider another func- 
tion of integrated level-spacing distribution defined by 
F(s) = J s °° P{s)ds = 1 - I P (s). The meaning of F(s) 
is the integrated fraction of level-spacings larger than s. 
Since F(s) of Pgue(s) is less than that of Ppe(s) at 
large s, we may expect that F(s) at larger s decreases 
or remains unchanged with increasing L for extended 
states while it increases with increasing L for localized 
states. FiglTTl shows curves of F(s) at (E = 0,P = 0.1) 
(a), (E = 0.02, P = 0.1) (b), and (E = 0.5, P = 1.5) 
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FIG. 11: F{s) vs. s for L = 8, 12, 16, 20, 24, and comparison 
with that of Pgue(s). (a) E = and P = 0.1; (b) E = 0.02 
and P = 0.1; (c) E = 0.5 and P = 1.5. 



FIG. 12: F{s) vs. s for L = 8, 12, 16, 20, 24 and comparison 
with that of Pque{s). (a) E = and P = 0.7; (b) £ = 0.02 
and P = 0.7; (c) E = 0.5 and P = 0.5. 



(c) for L — 8, 12, 16, 20, 24 and comparison with that of 
Pgue{s). FigH3 shows (E = 0, P = 0.7) (a), [E = 
0.02, P = 0.7) (b), and (E = 0.5, P = 0.5) (c). In view 
of FiglHl an d Fig^l it is clear that the results of F(s) 
coincide with those of Ip(s) concerning the localization 
property. We also calculate F(s = 2) for each fixed en- 
ergy in the whole range of mixing. As shown above, the 
same criteria as that for Iq and Ip(s = 0.5) can be em- 
ployed. The curves of F(s — 2) vs. P are plotted in 
Fig|T21for E = (a), E = 0.02 (b) and E = 0.5 (c). 
One can see that they are consistent with the results of 
I (Fig0 and I P (s = 0.5) (FigJTDJl. 



B. Discussion of the localization property 

In the last subsection, we analyzed the global shape 
of P(s) and its behavior at small and large s by con- 
sidering lo, Ip(s) and F(s), respectively. Analysis of all 
these quantities leads to essentially the same conclusion 
concerning the localization property, as follows. For zero 
interband mixing, only states at the two LB centers are 
extended. In the presence of interband mixing, new ex- 
tended states emerge. States near the LB centers, — i.e., 
E ~ 0, — are delocalized by weak interband mixing and 
localized by strong mixing, with a transition point at 
some intermediate mixing P c . For states near the region 
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FIG. 13: F(s = 2, P) vs. P for L = 8, 12, 16, 20, 24, (a) £ = 0; 
(b) E = 0.02; (c) E = 0.5. 



or P > P c = 1 for a square lattice^. Therefore, an ex- 
tended state is formed by strong mixing. One hopes that 
the intra-band tunnelings at nodes will only modify the 
threshold value of the mixing strength. 
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between two LBs, — i.e., E ~ 0.5, — they are localized at 
both weak and intermediate mixing and delocalized by 
strong mixing. 

The existence of new extended states at E ~ 0.5 in the 
case of strong interband mixing can be understood as fol- 
lows. Assume that the intra-band tunneling at nodes are 
negligibly weak for states of E ~ 0.5, we saw already from 
Fig. Ela) that the maximum interband mixing (sin# = 1) 
delocalizes the state, which is localized at zero interband 
mixing. If one views p = sin 2 9 as connection probability 
of two neighboring loops of opposite chirality, our two- 
channel model without intra-band tunnelings at nodes 
is analogous to a bond-percolation problem. It is well- 
known that a percolation cluster exists at p > p c = 1/2 



FIG. 14: J vs. E for L = 8,12,16,20, (a) P = 0.1; (b) 
P = 0.7; (c) P = 1.5. 

In order to show explicitly the existence of a narrow 
band of extended states and its evolution with increasing 
mixing, curves of Iq vs. E are plotted for three values 
of P in Fig^J A band of extended states is formed 
around the LB center E £ [0,0.1] for P = 0.1. When 
P is increased to an intermediate value 0.7, this band of 
extended states is lifted up to E £ [0.8, 1.6]. For strong 
mixing, it is further shifted to E G [0.4,0.5]. Thus the 
band of extended states in the lower LB emerges in weak 
mixing and tends to float up in energy with increased 
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mixing. By symmetry one can expect that the extended 
band of the upper LB should tend to dive down in energy 
with increased mixing. The two bands of extended states 
in the lower and upper LBs should finally meet at the 
middle energy region in the case of strong mixing. 

The above results are restricted to the case of two LBs, 
while there can be infinite LBs in the continuum model 
for a realistic system. Combine our above results and the 
fioat-up-merge picture proposed by Sheng et. aliiS, one 
can expect the following conclusions for multiple LBs. In 
the case of weak interband mixing, a narrow extended 
band emerges in each LB. With increasing mixing, i.e., 
increasing disorders or decreasing magnetic field, the ex- 
tended band in the lowest LB floats up and finally merges 
with that in the second lowest LB. Then, this extended 
band will further shift up and merge into that in the third 
lowest LB, and so on so forth. 

To express our numerical results in the plane of energy 
and interband mixing, a topological phase diagram shown 
in Fig. I15f a) is obtained. In the absence of interband 
mixing, only the singular energy level at each LB center 
is extended. In the presence of interband mixing of op- 
posite chirality, there are two regimes. At weak mixing, 
each of the extended states broadens into a narrow band 
of extended states near the LB centers. With increased 
mixing, the extended states in the lowest LB shift from 
the LB center(see Fig ll4|l . These extended states will 
eventually merge with those from the higher LBs. This 
shifting of extended states was also observed before^. At 
strong mixing, a band of extended states exists between 
neighboring LBs where all states are localized without 
the mixing. 

Let us look at the consequences of the above results. 
For weak disordered systems in IQHE regime, the Lan- 
dau gap is larger than the LB bandwidth. Thus there is 
no overlap between adjacent LBs. According to the semi- 
classical picture, electronic states between the two adja- 
cent LBs should be from either the upper or the lower 
bands with the same chirality in this case. It means that 
no interband mixing occurs and there is only one ex- 
tended state in each LB. This may explain why scaling 
behaviors were observed for plateau transitions in early 
experiments on clean samples. Interband mixings occur 
when the Landau gap is less than the LB bandwidth. Sys- 
tems of relatively strong disorders in IQHE regime should 
correspond to this case. As the single extended state at 
each LB center broadens into a narrow extended band, 
a narrow metallic phase emerges between two neighbor- 
ing IQHE phases. Thus each plateau transition contains 
two consecutive quantum phase transitions for strongly 
disordered systems. The bands of extended states will 
merge together in strong mixing. This strong mixing 
regime corresponds to the case when the Landau gap is 
much smaller than the bandwidth. Since the Landau 
gap is proportional to the magnetic field, the disordered 
system should always enter the strong mixing regime be- 
fore it reaches the weak field insulating phase, regardless 
of how weak the disorders are. In terms of QH plateau 



transitions, a direct transition occurs because a narrow 
metallic phase exists between two QH phases in a weak 
field. Thus, we propose that a direct transition from an 
IQHE phase to the insulating phase at weak field is real- 
ized by passing through a metallic phase, and it should 
hold for both weak and strong disordered systems. 
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FIG. 15: (a) Topological phase diagram of electron localiza- 
tion in E — P plane. The shadowed regime is for extended 
states (metallic phase), (b) Topological QH phase diagram in 
W ~ B plane. W stands for the disorder strength, and B for 
the magnetic field. The shadowed regime is for the metallic 
phase. The area indicated by the symbol n is the n-plateau 
IQHE phase. The rest area is the insulating phase. 

Plot above results in the plane of disorder and the mag- 
netic field, we obtain a new topological QH phase dia- 
gram as shown in Fig. I15f b). This is similar to the em- 
pirical diagram obtained experimentally in Ref. 14. The 
origin (W = 0, 5 = 0) is a singular point. According to 
the weak localization theorjii, no extended state exists 
at this point. Differing from existing theories, there ex- 
ists a narrow metallic phase between two adjacent IQHE 
phases and between an IQHE phase and an insulating 
phase. This new phase diagram is consistent with the 
non-scaling experiments 1 - where samples are relatively 
dirty, and interband mixing is strong, corresponding to 
a process along line a in Fig. I15f b). The system under- 
goes two quantum phase transitions each time it moves 
from the QH insulating phase to IQHE phase of n = 1 
and back to the weak field insulating phase as the mag- 
netic field decreases. To verify this claim, we analyzed 
the original experimental data in Ref. 17 according to 
the assumption of two quantum phase transition points. 
The experiment data of the logarithm of the longitudi- 
nal resistance ln[R xx (f,T)] are shown in Fig. 1161 where 
/ is the filling factor of LBs and T is the temperature. 
According to the theory of continuous transitions, one 
should obtain 

ln[Rxx{v,T)} = Fi(Si(f)/T) (12) 
with Si(v) — (fd - f) ZlVl for the region of / < f cl while 



ln[R xx (v,T)} = F 2 (S 2 (f)/T) 



(13) 



with Sa(y) ~ (/ - f C 2) Z2 " 2 for the region of / > f c2 . 
Previous theories predict one single critical point, — i.e., 
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FIG. 16: Experiment data of the logarithm of the longitudi- 
nal resistance \u\R xx {f,T)/R xx {Q.QA7,4QmK)\ in Ref. 17. / 
is the filling factor of LBs and T is the temperature. 



fci — fc2 and Z\V\ = z 2 v 2 . But our results suggest two 
distinct critical points. By standard scaling analysis, two 
good scaling behaviors are obtained for two close critical 
filling factors of f c \ — 0.6453 and f c2 = 0.6477 as shown 
in Fig. El The critical exponents in both the left side 
and the right side of the transition region are equal to 
the value zv — 2.33 ±0.01. On the other hand, the fit for 
one single critical point fails. Fig. shows the result of 
a single critical point at v c = 0.646. It is the best fitting 
result for a single critical point if we require that the two 
critical exponents are approximately equal and the scal- 
ing law is optimally obeyed. The two critical exponents 
are Z\v x = 2.58 ± 0.02 and z 2 v 2 = 2.60 ± 0.02, deviating 
from the theoretical results v ~ 2.33. One can also see 
clearly systematic deviations from the scaling law in the 
region close to the critical point at both sides in Fig. fT81 
This implies that the transition process is governed by 
two separated critical points instead of one. The regime 
between the two critical points should correspond to the 
metallic phase. Our fitting shows that the width of this 
regime is about 5 x I0~ 3 tesla while the value of the mag- 
netic field was increased by 1 ~ 2 x 10~ 3 tesla each time 
in the experiments. This may explain why the metallic 
phase was overlooked. 

It is worth noting that there is another puzzle in the 
non-scaling experiment which may be solved by our two- 
critical-point picture. To be specific, let us consider the 
experimental data for the transition between the QH 
insulating phase and the n = 1 IQHE phase. It was 
shown that the logarithm of the longitudinal resistance 
ln[R xx (f, T)] can be fitted with a linear function of the 
LB filling factor / (see Fig ll6|l as following^ 

ln[R xx (f,T)]=ln[R xx (f c ,T)]-(f-f c )/(a + 0T)] (14) 

where a and (3 are positive constants, f c is the filling fac- 
tor where curves of different temperature T cross approx- 
imately. Since a is non-zero^, it leads to the conclusion 
that R xx (f,T) at the limit of T = remains finite for 
every f. This is puzzling because it is inconsistent with 



FIG. 17: The fitting result of two critical points at the left and 
the right side. The two straight lines show coincidence with 
the scaling law. The critical filling factors are f c i ~ 0.6453 
and f C 2 = 0.6477. The two critical exponents are equal to the 
value zv = 2.33 ± 0.01. 
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FIG. 18: The best fitting result of one single critical point at 
the left and right side. The critical filling factor is f c = 0.646. 
The two straight lines illustrate systematic deviations from 
the scaling law at regions close to the critical point. The 
average values of the two critical exponents are ZiVi = 2.58 ± 
0.02 and 22^2 = 2.60 ± 0.02, respectively. 



the theoretical requirement that R XX (T = 0) = 00 in the 
QH insulating phase, i.e., / < f c , and R XX {T = 0) = 
in the n = 1 IQHE phase, i.e., / > f c . This puzzle may 
be solved as follows. Combine the linear relationship be- 
tween ln[R xx (f, T)] and / for fixed T with our picture of 
two critical points f c ± < f c2 , we expect 

ln[R xx {f,T)]=ln[R xx {U,T)}-{f-f cX )/(AiT* v ) (15) 

in the QH insulating phase, i.e., / < / c i, while 

ln[R xx {f,T)] = ln[R xx (f c2 ,T)]-(f-f c2 )/(A 2 T™) (16) 

in the n=l IQHE phase, i.e., / > f c2 , where A\ and A 2 
are positive constants, and z and v are critical exponents. 
It is clear that both R xx (f, T = 0) = 00 in / < f cl and 
Rxx{f,T — 0) = in / > f c2 are recovered. While a 
finite value of R xx (/, T = 0) in the region f c i < f < 
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f C 2 is consistent with our prediction of a metallic phase 
between the two critical points. 

Our model can also be used to describe spin-polarized 
systems. In this case, the two LBs are for spin up and 
spin down states. Indeed, two-channel CC models have 
been used before to simulate a spin-resolved problem 6 . 
In the presence of interband mixing, two distinct criti- 
cal points were obtained. They were related to the two 
extended states in the two subbands, which are shifted 
by the mixings. However, this study could not discern 
whether the states in between are extended or localized. 
Thus it was not clear whether the extended states in 
the presence of mixing are just the two points or form a 
band. In this sense, our results are consistent with those 
of early works. It is worth noting that the spin resolved 
problem is different from the non-scaling experiment. In 
the spin resolved problem, the energy region considered 
includes the centers of both the spin-up and spin-down 
subbands. Thus there are two near-degenerate extended 
states in this region in the absence of interband mixing. 
It is then natural to regard the two distinct critical points 
as two distinct extended states of the two subbands in the 
presence of mixing. However, in the non-scaling exper- 
iment, the region considered includes only the center of 
the lowest LB. Thus there is only one extended state in 
the absence of interband mixing. In this case, the two 
separate critical points may not be considered as two 
distinct extended states. It seems that the only suitable 
explanation is the existence of a band of extended states 
between the two points. 

One should also notice that two types of metallic 
phases have been studied extensively in the QH system. 
One is the composite Fermion state at the half-filling 
in the lowest Landau level (LL) and the other is the 
stripe state at the half-filled higher LLs. These states 
are formed by the Coulomb interaction effect in the high 
mobility samples. They are different from our metallic 
phase due to level mixing in the paper. Although we 
have not considered the electron-electron interactions in 
our study, there is no reason why the derealization ef- 
fect dof level mixing will be diminished by the Coulomb 
interaction. Of course, the interaction could change the 
level mixing effect's dependence on the magnetic field. 



CONCLUSIONS 



agram is consistent with non-scaling behaviors observed 
in recent experiments. A standard scaling analysis on 
experiment data in Ref. 17 supports our results. 
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FIG. 19: A two-channel network model of 2 x 2 nodes with 
periodic boundaries along both directions. Z's are the wave- 
function amplitudes on links. The notations are as follows. 
H and V stand for horizontal and vertical links, respectively. 
u (I) is for the upper (lower) LB. Si are SO(4) matrices de- 
scribing tunneling at nodes, and Mi are U(2) matrices for 
interband mixing. 



In conclusion, we find by numerical calculations within 
the network model that the single extended state at each 
LB center in the absence of interband mixing broadens 
into a narrow band of extended states when the effect 
of mixing of states of opposite chirality is taken into ac- 
count. With the decrease of magnetic field or increase 
of disorders, these extended bands further broaden and 
may merge together. Based on this, we propose a new 
phase diagram in which a narrow metallic phase exists 
between two neighboring IQHE phases and between an 
IQHE phase and an insulating phase. This new phase di- 



In this appendix, we explicitly construct the evolution 
matrix U for a 2 x 2 two-channel CC-network model as 
shown in Fig. ^| Periodical boundary conditions in both 
directions are imposed as explained in section III. Z-s are 
the wavefunction amplitudes on links. The notations are 
as follows. H and V stand for horizontal and vertical 
links, respectively, u (I) is for the upper (lower) LB. 
Si are SO (4) matrices defined in Eq. [3| describing the 
tunneling on nodes, and Mi are U(2) matrices defined in 
Eqs. [S] and describing interband mixing. From Fig. 1191 
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we can obtain 
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where 1 and 6 are the 2x2 identity and zero matrices, 
respectively. If we define 
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then the evolution equation is 
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The evolution operator U is 



Si; Hi = 
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where is the 4x4 zero matrix. It has the structure of 
Eq. El 
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